Noise-induced inhibitory suppression of malfunction neural oscillators 
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Motivated by the aim to find new medical strategies to suppress undesirable neural synchroniza- 
tion we study the control of oscillations in a system of inhibitory coupled noisy oscillators. Using 
dynamical properties of inhibition, we find regimes when the malfunction oscillations can be sup- 
pressed but the information signal of a certain frequency can be transmitted through the system. 
The mechanism of this phenomenon is a resonant interplay of noise and the transmission signal 
provided by certain value of inhibitory coupling. Analyzing a system of three or four oscillators 
representing neural clusters, we show that this suppression can be effectively controlled by coupling 
and noise amplitudes. 

PACS numbers: 05.45.-a, 05.40.Ca 



I. INTRODUCTION 

There are diseases in brain function (specially epilepsy) 
that involve a large population of neurons firing pulses 
synchronously in an oscillatory manner, in contrast with 
their typical excitable behavior. It is then important to 
understand the mechanisms that, taking advantage of the 
natural connectivity between neurons, might help in the 
suppression of these oscillations related to the malfunc- 
tion of the brain. It would also be desirable, besides 
this suppression, to restore the transfer functionality such 
that information could be transmitted through the neural 
network. 

Several solutions to prevent undesirable neural syn- 
chronization have been proposed. For instance, one could 
perform a permanent high-frequency stimulation P, 
To avoid disadvantages of this strategy (high energy con- 
sumption and brain adaptation) the concept of phase 
resetting j3| has been developed and extended to pop- 
ulations of interacting oscillators subjected to random 
forces as demand- controlled deep brain stimulation tech- 
niques 0,13 . A promising framework for the development 
of future medical methods has been suggested by several 
recent theoretical results such as delay feedback control of 
collective synchrony or noise-induced excitability 0. 
In the latter technique not only the resynchronization of 
oscillations is achieved but also a restoration of informa- 
tion transmission is possible. 

The aim of this paper is to present a new method that 
uses the interplay between inhibitory coupling and noise, 
two usual ingredients of living neural networks, in or- 
der to suppress undesirable oscillations but nevertheless 
being able to transmit an information signal. More pre- 
cisely, we investigate in detail the effects of the coupling, 
signal frequency and noise intensity on the oscillation 
suppression with signal transmission restoration in a lin- 
ear chain of three or four units which combines excitable 
and oscillatory neuron models. 



To model the neurons we use the archetype FitzHugh- 
Nagumo equations 01 which are a simple example of two- 
dimensional excitable dynamics and exhibit a Hopf bifur- 
cation from excitable to oscillatory behavior as a control 
parameter is varied. The dynamics of oscillatory and ex- 
citable systems near the Hopf bifurcation attracts large 
interest because in this region their sensitivity is great- 
est and the system is predestined for a reliable signal 
response or information exchange. The study of Coher- 
ence Resonance [3 [13 (also named as Stochastic Co- 
herence [U) and Stochastic Resonance in nonlinear ex- 
citable units 0, 0, Q, 01 arose a strong interest on 
this field. 

It is important to remark here that most of the previ- 
ous work was carried out considering FitzHugh-Nagumo 
units coupled through its activator variable (the one 
that represents the membrane potential), whereas brain 
neural networks may contain both activator and in- 
hibitory synaptic connections It is known that the 
inhibitory coupling between identical oscillators evokes 
many limit cycles of different periods and with differ- 
ent phase relations 0, which are stable in large re- 
gions of the control parameter space and is usually re- 
ferred to as "de-phasing" 0, |22(| or "phase-repulsive" 
|2ll | interaction. De-phasing was shown to be a source 
of multi-rhythmicity in different systems [2^ [25l | . 

With neural noisy elements, a de-phasing interaction of 
stochastic limit cycles (instead of deterministic ones) may 
result in the coexistence of spatiotemporal regimes selec- 
tively sensitive to external signal periods. In such sys- 
tems, noise plays at least two roles: first, it stimulates 
firing of stable elements and, thereby, their interaction 
during return excursions; second, it stimulates transi- 
tions between coupling-dependent attractors if the life- 
time thereof is sufficiently long. 

In this paper we extend our research on the influence 
of inhibitory coupling. In contrast to our previous in- 
vestigation of frequency selective Stochastic Resonance 
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in linear chains of identical excitable FitzHugh-Nagumo 
models |26i,|23|, we consider chains of nonidentical oscilla- 
tors and focus on the influence of the oscillatory units on 
the dynamics of the whole chain. Further on, we are in- 
terested in the suppression of the oscillations and in par- 
allel in an information exchange along the chain across 
the suppressed elements. 

The paper is organized as follows. In the next section 
II we start with an investigation of three inhibitory cou- 
pled units, each of them representing a cluster of neural 
oscillators. A similar architecture may be responsible, for 
example, for the activities of neural circuits in a nucleus 
found in the brains of songbirds [2^. In such circuits 
the connection between different functional units of the 
brain is mainly due to inhibitory coupling, whereas the 
connections within each unit are mainly through the acti- 
vator variables. Since strong activator coupling tends to 
synchronize the population of interacting units one can 
neglect, as a first approximation, that each functional 
unit is composed itself of several units, and restrict one- 
self to a case in which only one (mean) unit is considered 
for each region of the brain, coupled with others through 
the inhibitor variable. In this architecture the oscillat- 
ing element is directly surrounded by inhibitory coupled 
excitable elements. To extend the study onto another 
architecture, in section III, we analyze a circuit where 
two coupled oscillating units are connected from both 
sides with excitable element. Some general conclusions 
are drawn in the last section IV. 



II. 



THREE NON-IDENTICAL INHIBITORY 
COUPLED FITZHUGH-NAGUMO 
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FIG. 1: Scheme of the setup for the case A'^ = 3. At both ends 
there are excitable units, coupled through inhibitor coupling 
to the middle (oscillatory) element. 

We consider a rather simplistic model with a minimal 
scheme of connections that can retain the basic structure 
of the system we want to mimic: i.e. a signal injected 
through a neuron models the information to be trans- 
mitted along the network with one (or more) malfunc- 
tion oscillator. The basic setup is an open, linear chain 
where both ends have an excitable unit. The middle unit 
represents an ill (oscillatory) unit (Fig. QJ. We want to 
study whether a periodic, subthreshold signal, acting on 
the left element can reach the right one, in such a way 
that no oscillations appear in the middle unit. 



The scheme in Fig. ^ corresponds to three FitzHugh- 
Nagumo oscillators coupled through the inhibitory vari- 
ables: 



exi = 2/1 - y + (1) 

yi = ai - xi + £,i{t) + As cos (tut) + D{y2 - yi) 

£±2 = 2/2 - y + 2:2 (2) 

2/2 a2 - 2^2 + + -0(2/1 - 2/2) + £'(2/3 - 2/2) 



£2^3 = 2/3 - y + 2^3 
2/3 = as - ^3(t) + D(2/2 - 2/3) 



(3) 



Where lj — li^jTs, is the frequency of the input signal 
with period Tg . The Gaussian (white) noise sources (t) 

satisfy (6(0?j (i')) = ^^IK^ " 

In a neural context, yiit) represents the membrane po- 
tential of the neuron and Xi(t) is related to the time- 
dependent conductance of the potassium channels in the 
membrane "29^. The dynamics of the activator variable 
yi is much faster than that of the inhibitor indi- 
cated by the small time-scale ratio parameter e. It is 
well known that for \ai\ > 1 a single unit has a stable 
fixed point and presents excitable behavior: small per- 
turbations are followed by a smooth return to the fixed 
point, while a perturbation larger than a threshold value 
induces a return through a large excursion in phase space 
(a spike). For |aj| < 1, the fixed point becomes unstable 
and a stable limit cycle appears. In this regime, the dy- 
namics consists in a periodic series of spikes. Along this 
section, we will consider the fixed parameters: e = 10~^, 
01 = ^3 = 1-01 and ai = 0.99, such that the two end 
units are excitable and the middle one, oscillatory. We 
have checked that (in the absence of external forcing and 
noise) the three units retain their excitable or oscillatory 
character despite the coupling amongst them, such that 
the middle unit spikes periodically and the two end units 
display small subthreshold oscillations around the fixed 
point. 

The issue now is the behavior of these same units when 
noise and external forcing are present. We will show that 
it is possible to have a noise-induced regime in which 
the oscillations of the middle unit are suppressed. To 
characterize this phenomenon of oscillation suppression, 

we have computed ^3 -^s*\ the number of spikes per 
time unit at the i— th neuron, defined as the number of 

times the variable yi(t) surpasses a fixed threshold. Tvi'"* 
represents the reciprocal of the averaged inter-spike time 
interval. 

An important point is whether in this oscillation sup- 
pression regime, noise can help to transmit the informa- 
tion of the subthreshold external signal by a stochastic 
resonance mechanism. In order to address this issue, we 
compute the linear response, Q*-*-*, of the i— th neuron in 
the chain at the input frequency lo [sol l3ll | : 



Q« = (|22/,(t)e^"*d). 



(4) 
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where (...) denotes a time average. 



A. Oscillation suppression via a noise-induced 
dynamical trap 

In a previous work |32j| it has been shown that in a 
system of three FitzHugh-Nagumo units in the oscilla- 
tory regime (and in the absence of external forcing) the 
inhibitory coupling leads to two coexisting dynamical at- 
tractors, with different natural frequencies. These attrac- 
tors correspond to an anti-phase oscillator movement and 
to the so-called dynamical trap regime where the middle 
oscillator is at rest and the two oscillators at the ends 
oscillate in anti-phase. If one now applies a weak ex- 
ternal periodic signal to one of the end units, one can 
still achieve the suppression of the self-excited oscilla- 
tions of the middle unit and, at the same time, achieve 
a reliable transmission of the signal, provided the follow- 
ing two conditions hold: (i) the frequency of the external 
signal coincides with the natural frequency of the dynam- 
ical trap attractor, and (ii) the noise in the system is near 
the optimal one for the desired signal amplification (i.e. 
stochastic resonance phenomenon on this attractor). 

A similar result appears in our system of three coupled 
units. For very small noise, the situation is as described 
at the beginning of the section with the middle unit os- 
cillating and the end units at rest. As noise increases, 
one observes random switches between this state and a 
dynamical trap regime in which the middle unit is at rest 
and the two end units spike in anti-phase. This effect can 
be quantified by measuring the number of spikes ivi*'' and 
the responses Q*^*' as a function of the noise intensity. As 
shown in Figs. El one can distinguish several behaviors 
depending on the period of the external forcing. 

(a) This is the case where the period of the input sig- 
nal equals the natural period of an isolated FitzHugh- 
Nagumo oscillator (T, = 2.8 for a — 0.99). The noise- 
induced oscillation suppression described before is appar- 
ent in the right panel of this figure, where it is shown that 

(2) 

the number of spikes at the middle unit, Ns , first de- 
creases as the noise intensity increases. This oscillation 
suppression is maximum at a value of the noise intensity, 
CT^ « 2-10^^. At noises larger that this value, the number 
of spikes in the three units are very close to each other. 

In the left panel we plot the response Q^'^ of each unit. 
Note that there is a range of values for the noise inten- 
sity for which the middle unit responds to the injected 
signal most effectively than the end units, as signaled by 
a higher value of the response Q^^^. For increasing noise 
intensity, beyond the value where the oscillation suppres- 
sion was maximum, all units have a similar response. 

(b) For an intermediate range of periods Tg G [3,3.4], 
we observe that there exists a range of noise intensities 
(cr^ G [10~^,10~^]) such that the number of spikes is 



strongly reduced in the middle oscillatory unit, while the 
response to the driving frequency is better than in the os- 
cillatory unit, i.e. this is the manifestation of the dynamic 
trap regime. One can clearly see the effective oscillation 
suppression of the oscillatory middle element (the mal- 
fmiction neuron, see Fig. right) and -despite of this 
suppression in the middle of the chain- the reliable infor- 
mation transport from the first to the last neuron model 
by a large linear response Q in these elements (Fig. EJ), 
left). The dynamic trap regime includes an anti-phase 
motion of the first and the last unit which results in com- 
bination with the inhibitory coupling in a suppression of 
the oscillations of the middle element. 

(c) Increasing even further the period, Tg — 4.5, the 
external signal is now in resonance with, and amplifies, 
the anti-phase motion in which the first and the last 
units oscillate in-phasc and in anti-phase with the mid- 
dle one. In this case, another interesting regime appears 
in the noise range cr^ G [10^^, 10^^] as observed in the 
spike number diagram (right panel of Fig. 13;), where the 
spike numbers of all three elements coincide nearly, as 
well as in the linear response plot (left panel), where 
all oscillators display a very similar linear response Q. 
This anti-phase regime demonstrates a totally different 
behavior than the dynamic trap regime, case (b) dis- 
cussed previously. Note that the anti-phase regime ap- 
pears for a much larger noise intensity than the dynamic 
trap regime, hence showing a double selectivity by the 
input frequency and the noise intensity. 

(d) Finally, for much larger period, = 6.0, there is no 
resonance, Figs. This can be observed especially at 
the linear response Q (left panel) which is much smaller 
than in the resonant cases. Noteworthy, the last element 
in the chain exhibits a poor signal response. 

B. Control of suppression by the coupling strength 

Noise-induced dynamical trap suppression is possible 
by the existence of a new attractor originated in the in- 
hibitory nature of the coupling. Hence, the coupling in- 
tensity, D, controls the effectivity of the suppression, as 
well as the frequency of the attractor. The existence of 
an optimal value for D is shown in Fig. |31 in which we 
plot the linear response Q and the spike numbers Ns as 
a function of the coupling in intensity. Starting from the 
maximum of the resonance curve of the linear response 
Q on the noise (Fig. I^b right panel) , it is clearly observed 
that there is an optimal value of the inhibitory coupling 
D such that the middle unit is silent (Fig. Alright panel), 
while the first and last units effectively respond to the 
driving frequency (Fig. |31 left panel). 

Since both types of coupling, inhibitory and activator, 
can be immanent in neural networks, we have investi- 
gated how the suppression can be regulated if we tune 
the coupling from an activator to an inhibitory one. To 
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FIG. 2: Optimal noise suppresses malfunction oscillations 
while letting the signal (within a certain range) to be trans- 
mitted. This effect occurs due to dynamical trap, supported 
by inhibitory coupling, (a) non-resonant, Ts = 2.8; (b) dy- 
namic trap, Ts — 3.1; (c) anti-phase resonance, Ts — 4.5; 
(d) no resonance, Ts = 6.0. Other parameters: e = 10~^, 
ai,3 = 1.01, 02 = 0.99, As = 0.01, D = 0.15. The left and 
right columns correspond to the Q and Ns measures. 
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FIG. 3: Linear response Q (left), and normalized spike num- 
ber Ns (right), versus inhibitory coupling. The other param- 
eters are: e = 10"", m.a = 1.01, a2 = 0.99, As = 0.01, 



Ts=3.1 and a„ 



3 • 10" 



inhibitory coupling (from right to left) leads an abrupt 
suppression of the middle oscillator (Fig. 0] left) and to 
a significant joint increase of the linear response Q of 
the first and third oscillators, but not of the middle one 
(Fig. fright). Note the logarithmic scaling of the param- 
eter a at the abscissa. We clearly observe, that already a 
small fraction of activator coupling (in the order of 1%) 
destroys the dynamic trap regime in the given parameter 
set. 
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FIG. 4: The influence of the type of coupling on the linear 
response Q (left) and the normalized spike number Ns (right) 
for three coupled FitzHugh-Nagumo's. The sliding parameter 
a shifts the weight of the diffusion constant D from a pure 
inhibitory coupling (a — 0.0) to a pure activator coupling 
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0.99, < = 4 ■ 10-^ As = 0.01, r, = 3.1 and D = 0.15. 



do this, we have added activator coupling in the model by 
interchanging D hy {1 — a)D in the equations for the in- 
hibitory variable yi and inserting the terms aD {x2 — xi) 
in Eq.^ a D (xi—X2)+aD{x3~X2) in Eq.[21 aZ)(x2 — 2:3) 
in Eq. |31 These extensions of the model are used only in 
this section for the calculation of Fig. ^ With help of the 
new sliding parameter a we change the weight of the type 
coupling from a = I (pure activator coupling) to a = 
(pure inhibitory coupling). The results are illustrated in 
Fig. ^ We clearly see that increasing the weight of the 



III. MODEL FOR FOUR NON-IDENTICAL 
UNITS 

Next we delve into the question of whether larger 
chains with more coupled oscillatory units also show the 
same phenomenon discussed in the previous section. Al- 
though it would seem a rather trivial proposal just en- 
large it to a case in which the system size is iV = 4, 
the dynamical regimes that arise in such situation are 
far from being simple modifications of the results shown 
above. 
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We will not consider an enlargement of the excitable 
ends of the chain, because it is a well known fact that cou- 
pling them through the activator variable with a strong 
enough bind, will result in an entrainment of such sub- 
chain and the dynamical evolution of such units will be 
effectively that of one oscillator. Then, the most inter- 
esting question arises from the enlargement of the mid- 
dle part, that is composed by oscillatory (malfunctioning) 
units. So, we add to the scheme of three elements (Fig.QJ 
an oscillatory element in the middle position and couple 
it by an activator coupling with the other (identical) os- 
cillatory element and with an inhibitory coupling to the 
adjacent excitable element (Fig. [SI). 
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FIG. 5: Scheme of the setup for the case A'' = 4. While 
at both ends there are excitable units, in the center there 
are oscillatory ones. The coupling between units of different 
nature is inhibitory, and the coupling between the oscillatory 
units is through the activator variable. 

The mathematical description of the scheme in Fig.Elis 
given in Eas. l5l8l The two oscillatory units (malfunction 
neurons) are placed at the middle position and are both 
coupled to their adjacent excitable one by an inhibitory 
coupling as in the chain of three elements, whereas an 
activator coupling is set between them. As in the pre- 
vious section, independent additive white noises act on 
the units and an external, subthreshold, periodic signal 
drives only the first element: 
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We will fix along the following simulations the param- 
eters: 02,3 = 0.99 (oscillatory regime), 01^4 = 1.01 (ex- 
citable regime) and the signal intensity Ag = 0.01 (sub- 
threshold) . 

We are interested in the signal penetration along the 
chain from the first to the last element as a function 



of the signal period and the noise intensity. In order 
to investigate whether the same phenomenon appears in 
this chain, two different cases are considered: first, we 
take the optimal parameters from the case A'^ = 3 and 
make the coupling between the oscillatory units strong 
enough such that they become entrained. In the second 
case, we use a weaker activator coupling. 



A. Strong inter-oscillatory coupling 

Let us focus first on a regime of strong coupling among 
the oscillatory units. We use the following set of param- 
eters £ = 10"'', ai,4 = 1.01, 02,3 = 0.99, C = 0.80 and 
D = 0.22. In this case, and without an external peri- 
odic signal {As = 0.0) injected nor noise (cr^ = 0.0), the 
analysis of the power spectrum exhibits that the natu- 
ral period of the system is Tnat ~ 2.67. The oscillatory 
units exhibit their periodic oscillations at their natural 
frequency. The excitable units, at their time, show only 
subthreshold oscillations at the natural frequency of the 
oscillatory units. 

In the presence of an external signal Fig. illustrates 
the normalized spike number and the linear response Q 
as a function of the noise intensity for different driving 
periods Tg. Fig.|B^ depicts the results when the system is 
subjected both to noise and external signal and the signal 
period Tg = 2.61 is slightly below the natural period. It is 
observed (as in the = 3 case) that now the oscillatory 
units respond but not the excitable ones. 

Increasing Ts well over the natural frequency, e.g. 
Ts — 2.8 or 2.9 (Figs.Eb orEt), the dynamic trap regime 
appears. It is important to emphasize that the quality of 
the signal transmission to the last unit is enhanced with 
respect to the = 3 case (compare Figs. 12 and Figs.lHll. 

An interesting phenomenon occurs for = 2.9 
(Fig. Et) where there are two well-differenced situations 
of dynamics trap like regimes. First, for very low noise 
intensities (cr^ « 10"^) there is an almost perfect sup- 
pression of the oscillations and at the same time a per- 
fect signal transmission which is the result of the desired 
dynamic trap regime. There is a then a secondary os- 
cillation suppression regime at cr^ « 2.5 • 10^^, at which 
the signal is not transmitted with the same fidelity as 
compared to the case at about « 10~^. What is hap- 
pening in the second regime is that the last unit is oscil- 
lating, neither at the driving frequency, nor at the natural 
one of the middle oscillators, but at another one. Fig. 
shows the power spectrum for such secondary regime in 
the interesting frequency range around driving and res- 
onance frequency. Let us consider the particular case of 
the fourth oscillator. It is subject to two different signals, 
one of them with the natural frequency of the third unit, 
and one with the external driving frequency. It is not 
trivial how this two signals interact in order to produce 
this unit's response, but it has been demonstrated that 
in non-linear systems [s^l subjected to two signals, the 
response may appear at neither any of the driving ones. 
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Nevertheless these facts, the important footprint of this 
secondary regime is the low response of the last unit to 
the driving frequency. 

One can clearly see in Fig. |7| three peaks in the fre- 
quency range uj G [1.8,2.8] in the system output. The 
first and highest peak at a; « 2.16 is well pronounced 
only for the first and driven oscillator and corresponds 
to a period T « 2.9, equal to the driving period Tg, i.e. 
only the driven oscillator exhibits a good response to the 
signal. The second peak, very close to the first one, at 
UJ « 2.245 (T « 2.8) is displayed mainly by the last unit. 
The third peak at a; 2.49 (T « 2.52) can be found in aU 
elements with nearly equal hight. The third peak has a 
very small influence on the total responses of the system 
(note the logarithmic scale) and it is produced by small 
sub-threshold oscillations. The corresponding time series 
is as follows: the first (driven) oscillator shows a reliable 
spiking behavior with a period equal to the driving pe- 
riod, while the two middle oscillators are mostly silent 
and the last oscillator spikes with a slightly reduced pe- 
riod of Tg = 2.8 (so leading to the difference in the linear 
response at the signal frequency Q between the first and 
last FitzHugh-Nagumo in this regime) . Due to the small 
difference in the periodicity, there is no phase locking in 
this regime and a continuous phase slip between the first 
and last unit appears. If the phase difference is large 
enough, the chain switches to an anti-phase regime, i.e. 
the otherwise silent middle oscillators spike in anti-phase 
to their excitable neighbors. This transition to the anti- 
phase attractor induces a delay of the last unit compared 
to the first one. This anti-phase regime is unstable at the 
considered parameter set and the chain switches back to 
the previous attractor with the silent middle elements 
and the phase slip between the first and last one. The 
interruption of this long-life attractor by the unstable 
anti-phase attractor results in a nearly equal spike num- 
ber of the first and last unit. Therefore, the interesting 
behavior in Figs.|Ht at a noise intensity cr^ « 2.5 • 10^^ is 
caused by a regime which is only similar to the dynamic 
trap regime, but is not exactly the desired dynamic trap 
and hence does not provide a reliable information trans- 
port. 

As a summary of this section, it could be said that 
the dynamic trap regime still occurs, but in a narrower 
region of the driving period (Ts G [2.8,3.0], Fig.lHJ than 
in the case of iV = 3 (T^ £ [3.0, 3.4]). 

Finally, Fig. |51 yields that there is also a range of 
inhibitory coupling D such that this phenomenon hap- 
pens. This resonance like behavior with respect to the 
inhibitory coupling strength is caused by the influence 
of this parameter on the resonance frequency of the dy- 
namic trap regime. This figure shows the existence of a 
maximum (located at a coupling D « 0.25) in the re- 
sponse as a function of this parameter. 




FIG. 6: Linear response Q (left column) and Ns (normalized 
spike number), in the right column, versus noise intensity for 
a chain of four oscillators, (a) Ts = 2.61; (b) Ts — 2.8; (c) 
Ts = 2.9; (d) Ts = 3.1; Other parameters: e = 10"*, aiA = 
1.01, 02,3 = 0.99, As = 0.01. The couplings are, C = 0.80 
and D = 0.22 (strong coupling regime), and Ts = 2.61. 



B. Intermediate inter-oscillatory coupling 

We also found another kind of dynamic regime in this 
model with a smaller activator coupling. The analysis of 
the power spectrum in the absence of noise and external 
signal shows that the oscillatory units exhibit their peri- 
odic oscillations at their natural period Tnat ~ 2.54. The 
excitable units, at their time, generate only subthreshold 
oscillations at the natural frequency of the oscillatory 
units. Note that the natural frequency is shifted from 
the previous case of a strong inter-oscillatory coupling 
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m 

FIG. 7: Power spectrum of a four oscillator system for the 
parameters: e — 10^*, ai.4 = 1.01, 02,3 = 0.99, As = 0.01, 
Ts = 2.9, C = 0.80, D = 0^22, and cr^ = 2.56 ■ 10"^ 
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FIG. 8: Linear response Q and normalized spike number Ns 
versus time periodicity Tg. The system is composed by four 
units in the strong-coupling regime and the rest of parameters, 
are: e = 10"*, 01,4 = 1.01, aa.a = 0.99, = 0.01, C = 0.80, 
D = 0.22, andaa = 2-10"^ ' 




D D 

FIG. 9: Linear response Q and normalized spike number Ns as 
a function of the inhibitory coupling strength D. The system 
is composed by four units in the strong-coupling regime and 
the rest of parameters, are: e = 10~*, ai.4 — 1.01, a2,3 = 0.99, 
As = 0.01, C = 0.80, cr^ = 2 ■ 10^^ and Ts = 2.9. 



{Tnat ~ 2.67). In this case, however, the dependence 
with Ts of the linear response Q curves and the oscilla- 
tion suppression is quite different from the previous case. 
Even for sUghtly detuned input signals Ts — 2.55 a strong 
dynamic trap arises in the system fFig. llOb). This oscil- 
lation suppression mechanism is very robust over a wide 
range of driving periods Ts fFigs. ITUk-c'). whereas a reli- 
able signal transmission along the chain can be observed 
only in a much narrower range of the driving period, 
Ts e [2.6,2.65] (Fig.Cnb). 

The oscillation suppression here is really robust, show- 



ing that the middle units do not spike for very large pe- 
riods of time. The phenomenon is also robust to changes 
of almost four decades in the noise intensity. 

Note that the curves for the number of spikes show an 
exact coincidence between the first and the last units (i.e. 
the excitable ones), although such a perfect matching 
does not occur for the linear response Q. The first and 
the last units fire at the same rate (same normalized spike 
number Ns), but they are not phase locked, i.e. there 
is a random phase slip. When the difference in phase 
between these two excitable units is large enough, this 
dynamic regime destabilizes and a regime in which there 
is an in-phase motion of excitable units, and (in anti- 
phase) spikes of the oscillatory units appears. But this 
last dynamic regime is unstable and rapidly falls to the 
previous one. It is interesting that the matching in the 
number of spikes occurs in the dynamic trap regime, i.e. 
that the sub-threshold dynamics of the oscillatory units is 
sufhcient to carry information from one end of the chain 
to the other one. 




FIG. 10: Linear response Q and normalized spike number Ns 
(left and right columns, respectively) versus noise intensity. 
The time periodicities, are (a) Ts = 2.55; (b) Ts = 2.61; (c) 
Ts — 5.2. The other parameters are: e = 10~*, 01,4 = 1.01, 
a2,3 = 0.99, As = 0.01, C = 0.20, and D = 0.50. 

Figs.^Jand^ldemonstrate that there are optimal val- 
ues of couplings for the suppression to occur. While the 
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dependence on the activator coupling C is such that the 
suppression holds for couplings larger than a given value, 
we observe a much narrower range, a resonance-like be- 
havior, as a function of the inhibitor coupling D. Even 
further, for D large enough, the oscillation suppression 
phenomenon disappears, and most of the spikes occur at 
frequencies different to the driving one (i.e. Q vanishes). 

Fig. 1131 shows the dependence on the signal periodic- 
ity Ts- It is clear that the oscillation suppression and 
signal transmission are optimal at the same parameter 
values. Furthermore, in the same figure it can be seen 
that there is a very narrow peak around the natural pe- 
riod {Ts = Tnat = 2.54) of the oscillatory units at which 
they respond optimally. Note that the oscillation sup- 
pression holds for a wide range of values of the driving 
period Ts- But the main result shown in this figure is the 
fact that the suppression of oscillations in the oscillatory 
units is much more robust than in the previous cases, i.e. 
TV = 3 and = 4 with strong coupling. This result is 
somewhat unexpected given the fact that the couplings 
are not as strong as in the previous parameter sets, and 
then the units are allowed to move more freely. 



(a) 0.5^ — — ^ — ^ — ^ — -| ' ' — ~ — — ^ — ^ — ^0.4 




0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 

c c 



FIG. 11: Linear response Q and normalized spike number Ns 
versus activator coupling C. The system is composed by four 
units, and the other parameters are: e — 10"*, ai,4 = 1.01, 
a2,3 = 0.99, A, = 0.01, D = 0.50, = 2.61 and cr^ = 2-10"^ 



(b)0.5i I i0.4 




D D 

FIG. 12: Linear response Q and normalized spike number Ns 
versus inhibitory coupling D. The parameters are: e — 10"'*, 
ai,4 = 1.01, a2,3 = 0.99, As = 0.01, C = 0.20, Ts = 2.61 and 

To show the different influence of activator and in- 
hibitory couplings, we have added an extra activator cou- 
pling in the model by interchanging D by {1 ~ a)D in 
the equations for the inhibitory variable yi and inserting 



•-.,■=1 









2.6 2.7 2.8 2.9 3.0 2.6 2.7 2.8 2.9 3.0 



T T 

FIG. 13: Linear response Q and normalized spike number 
Ns versus driving period Ts, For this four units system, the 
parameters are: e = 10~^, ai,4 = 1.01, 02,3 — 0.99, As — 0.01, 
C = 0.20, D = 0.50 and al = 2- 10"^ 

the terms aD{x2 — xi) in Eq. [Sj aD{xi — X2) in Eq. (HJ 
aD{x4 — X3) in Eq. [7| aD{x^ — X4) in Eq. |S1 These ex- 
tensions of the model are used only in this section for the 
calculation of Fig. ^1 We shift the balance between the 
activator and inhibitory coupling between these elements 
continuously with the parameter a. As shown in Fig. 1141 
a very sharp transition to a situation of oscillation sup- 
pression and no response to the driving frequency in the 
middle units is observed when the activator coupling is 
strong enough, a > 3 • 10"'^). Fig. E| (as Fig. ^for the 
= 3 case) shows the essential imperative of the in- 
hibitory coupling between the excitable and oscillatory 
units to reach the dynamic trap regime with the desired 
feature of oscillation suppression and information trans- 
mission. 




FIG. 14: The influence of the type of coupling on the the 
linear response Q and the normalized spike number Ns in a 
system composed by four coupled FitzHugh-Nagumo's. The 
sliding parameter a shifts the weight of the difi'usion constant 
D from a pure inhibitory coupling (a — 0.0) to a pure acti- 
vator coupling (a = 1). The other parameters are e = 10"'*, 
ai,4 = 1.01, a2,3 = 0.99, al = 2 ■ 10"^ As = 0.01, Ts = 2.61, 
C = 0.20 and D = 0.50. 



IV. CONCLUSIONS 

By studying coupled FitzHugh-Nagumo's units (or 
clusters, when one FitzHugh-Nagumo model represents 
many similar neurons), we have found new mechanisms 
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to suppress undesirable self-sustained oscillations in neu- 
ron networks with a reconstruction of the desired func- 
tionality to submit information with the help of a con- 
structive role of noise. We have shown that an inhibitory 
coupling between oscillatory (malfunction) and excitable 
neurons is essential to reach the dynamic trap regime 
which is responsible for the oscillation suppression and 
the information transport. The dynamic trap regime is 
characterized by an anti-phase spiking behavior of the 
excitable neurons at the end of the chain with the signal 
frequency and a silent (oscillation suppressed) behavior 
of the originally oscillatory units in between. The de- 
sired dynamic trap regime is sensitive with respect to 
the driving frequency, the noise intensity and the cou- 
pling strength. 

Besides the dynamic trap regime we found other at- 
tractors which are similar to the dynamic trap attrac- 
tor and offer also a reliable oscillation suppression but 
do not provide a good information transport along the 
chain. The oscillation suppression itself could be found 
in a larger parameter range than the combination of os- 
cillation suppression and reliable information transport. 

It is interesting to note that the suppression is present 
not only due to noise but also in presence of a strong 
enough driving force. It would then be interesting to 



study the effect of an alternative source of suppression of 
oscillations by the injection of an external signal. 

In the present work we have considered chains of three 
or four units to keep the number of parameters small. 
Our model could be exemplary also for larger systems, if 
one regards one oscillator in the model as a representation 
of a cluster of many neurons in a close area with similar 
properties. 

We have considered only paradigmatic models in a very 
general framework, hence we expect that these results are 
also relevant to other models with inhibitory coupling, 
which are used, for example, to describe various physical 
[s^ l- electronic chemical systems, biological 

systems with spatial non uniformities |37j . animal coat 
pattern formation js^, or artificial gene networks syn- 
chronization with slow auto- inducer diffusion j39l |40| . 
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